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Abstract. 

For finite quantum many-particle systems, a given system, induced by a transition operator, 
makes transitions from its states to the states of the same system or to those of another system. 
Examples are electromagnetic transitions (then the initial and final systems are same), nuclear 
beta and double beta decay (then the initial and final systems are different), particle addition 
to or removal from a given system and so on. Working towards developing a complete statistical 
theory for transition strength densities (transition strengths multiplied by the density of states 
at the initial and final energies), we have started a program to derive formulas for the lower 
order bivariate moments of the strength densities generated by a variety of transition operators. 
In this paper results are presented for a transition operator that removes ko number of particle 
by considering m spinless fermions in N single particle states. The Hamiltonian that is k- 
body is represented by EGUE(fc) [embedded Gaussian unitary ensemble of fc-body interactions] 
and similarly the transition operator by an appropriate independent EGUE. Employing the 
embedding U(N) algebra, finite-A formulas for moments up to order four are derived and 
they show that in general the smoothed (with respect to energy) bivariate transition strength 
densities take bivariate Gaussian form. Extension of these results to particle addition operator 
and beta decay type operators are discussed. 


1. Introduction 

Let us begin with the statement from the preface to the proceedings of the meeting held 
in 2006 on ’’Applications of Random Matrices in Physics”: [1] Random matrices are widely 
and successfully used in physics for almost 60-70 years, beginning with the works of Wigner 
and Dyson. Initially proposed to describe statistics of levels in complex nuclei, the Random 
Matrix Theory has grown far beyond nuclear physics, and also far beyond just level statistics. 
It is constantly developing into new areas of physics and mathematics, and now constitutes 
a part of the general culture and curriculum of a theoretical physicist. Besides applications 
in all branches of quantum physics, RMT is being used in disciplines such as Econophysics, 
Wireless communication, information theory, multivariate statistics, number theory, neural and 
biological networks and so on. The focus in this article is on the frontier topic of statistical 
properties of isolated finite many-particle quantum systems [2]. Examples for these systems are 
atoms, atomic nuclei, mesoscopic systems (quantum dots, small metallic grains), interacting spin 
systems modeling quantum computing core, ultra-cold atoms and so on. A route to investigate 
statistical properties is to employ the classical GOE or GUE or GSE random matrix ensembles 
with various deformations. For these, as Wigner states: the assumption is that the Hamiltonian 


which governs the behavior of a complicated system is a random real symmetric or complex 
Hermitian or Quaternion real matrix, with no special properties except for its symmetric or 
Hermitian or Quaternion real nature. However, for most of the isolated finite many-particle 
quantum systems, their constituents predominantly interact via two-particle interactions and 
the classical random matrix ensembles are too unspecific to account for this most important 
feature. One refinement which retains the basic stochastic approach but allows for this feature 
consists in the use of embedded random matrix ensembles [3, 4, 5, 6, 7, 8]. 

Therefore, it is more appropriate to represent an isolated finite interacting quantum system, 
say with m particles (fermions or bosons) in N single particle (sp) states by random matrix 
models generated by random A;-body (note that k < m and most often we have k = 2) 
interactions and propagate the information in the interaction to many particle spaces. Thus 
we have random interaction matrix models for m-particle systems. In the simplest version, 
the /c-particle Hamiltonian (H) of a spinless fermion (or boson) system is represented by 
GOE/GUE/GSE and then the m particle H matrix is generated using the m-particle Hilbert 
space geometry. The key element here is the recognition that the U(N) Lie algebra transports the 
information in the two-particle spaces to many-particle spaces. As a GOE/GUE/GSE random 
matrix ensemble in two-particle spaces is embedded in the m-particle H matrix, these ensembles 
are generically called embedded ensembles (EE). With GUE embedding we have EGUE and in 
this paper EGUE is used throughout. For EE, a general formulation for deriving analytical 
results is to use the Wigner-Racah algebra of the embedding Lie algebra [8]. The focus in the 
present paper is on transition strengths and they are not yet studied in any detail using EE [8]. 
Also, as emphasized in [4, 9], there are many open questions in the random matrix theory for 
transition strengths in finite interacting quantum many-particle systems. 

For finite quantum many-particle systems, induced by a transition operator, a given system 
makes transitions from its states to the states of the same system or to the states of another 
system. Examples are electromagnetic transitions (then the initial and final systems are same), 
nuclear beta and double beta decay (then the initial and final systems are different), particle 
addition to or removal from a given system and so on. Given a state with energy E t and say 
it is connected to a state Eg by a transition operator O, then the transition strengths are 
|(Ef | O | Ej)\ 2 and these will determine for example the life time of a state. It is important to 
recognize that the transition strengths probe the structure of the eigenfunctions of a quantum 
many-body system and thus they are very important from the point of view of experiments 
probing the structure of a system. Also, they are needed in many applications (for example, 
beta decay transition strengths are essential for nucleosynthesis studies). In the statistical 
theories, it is more useful to deal with the corresponding transition strength density (this will 
take into account degeneracies in the eigenvalues) defined by 

Io(Ei,E f ) = I(E f ) \{E f | O | Ef )| 2 I(Ef) . (1) 

In Eq. (1), 1(E) are state densities normalized to the dimension of the m particle spaces. Note 
that Ei and Ef belong to the same m particle system or different systems depending on the 
nature of O, the transition operator. 

Working towards developing a complete statistical theory for transition strength densities 
(transition strengths multiplied by the density of states at the initial and final energies) for 
isolated finite many-particle quantum systems, we have started a program to derive formulas for 
the lower order bivariate moments of the strength densities generated by a variety of transition 
operators. In general, the Hamiltonian may have many symmetries with the fermions (or bosons) 
carrying other degrees of freedom such as spin, orbital angular momentum, isospin and so on. 
Also we may have in the system different types of fermions (or bosons) and for example in atomic 
nuclei we have protons and neutrons. In addition, a transition operator may preserve particle 
number and other quantum numbers or it may change them. Among all these various situations, 



in our program we have considered five different systems: (i) a system of m spinless fermions 
and a transition operator that preserves the particle number; (ii) a system of to spinless fermions 
and a transition operator that removes say ko number of particles from the to fermion system; 
(iii) same as (ii) but for particles addition operator; (iv) a system with two types of spinless 
fermions with the transition operator changing ko number of particles of one type to ko number 
of other particles as in nuclear beta and double beta decay; (v) same as (i)-(iii) but for spinless 
boson systems. In [10], results are presented for (i). In the present paper results are presented 
for (ii) and their extensions to (iii) and (iv) are briefly discussed; Ref. [11] gives results from a 
first study for (iv). Let us add that (ii) and (iii) are important for example in nuclear physics 
as one and two particle removal and addition to a nucleus are important experimental probes of 
the structure of the atomic nucleus; see for example [12, 13]. Thus, the results in Sections 3-5 
have applications in nuclear physics. Now we will give a preview. 

Section 2 gives some basic results for EGUE(/c) for spinless fermion systems as derived in [14]. 
Using these results, formulas for the lower order bivariate moments of the transition strength 
densities for the situation (ii) above are derived and they are presented in detail in Section 3. 
Using these, results in the asymptotic limit are derived and they are presented in Section 4. 
Extension of the results in Section 3 to the situations (iii) and (iv) above are briefly discussed 
in Section 5. Finally, in Section 6 gives conclusions and future outlook. 

2. Basic EGUE(fc) results for a spinless fermion system 

Let us consider to spinless fermions in N degenerate sp states with the Hamiltonian H a &-body 
operator, 

H = Y f V ij {k) A\{k)Aj(k) , Vij{k) = (k,i\H\ k,j ) . (2) 

TO 

Here A'-(k) is a A: particle (normalized) creation operator and Ai(k ) is the corresponding 
annihilation operator (a hermitian conjugate). Also, i and j are fc-particle indices. The k 
and to particle space dimensions are (j() and (^) respectively. Representing the V matrix, 
defined by the matrix elements V t] , by GUE we have EGUE(/c) for the H matrix in m-particle 
spaces. Note that, for V a GUE, the real and imaginary parts of V t] are independent zero 
centered Gaussian random variables with variance satisfying, 

V ab (k) VcdJkj = 5 ad 6 bc . (3) 

Here the ’over-line’ indicates ensemble average. From now on we will drop the hat over H and 
denote when needed H by H(k). In physical systems in general k = 2 but in some systems such 
as atomic nuclei and BEC it is possible to have k = 3 and even k = 4 [15, 16, 17]. 

The U(N) algebra that generates the embedding, as shown in [14], gives formulas for the 
lower order moments of the one-point function, the eigenvalue density 1(E) = (( 6(H — E ))) and 
also for the two-point function in the eigenvalues. Used here is the U(N) tensorial decomposition 
of the H(k) operator giving v = 0,1,..., k irreducible parts B l,,Uv (k) and then, 

k 

H(k)= ]T W^(k)B^(k). (4) 

V=0\U}v£l/ 

Note that u v are labels of the irreducible representations (irreps) of the subalgebras of U(N) 
and their explicit structure will not play any role in the present work. With the GUE(A) 
representation for the H(k ) operator, the expansion coefficients W’s will be independent zero 
centered Gaussian random variables with 

W^(k)W U2 ^ 2 (k) = Vl 6^6^ 


( 5 ) 



For deriving formulas for the various moments, the first step is to apply the Wigner-Eckart 
theorem for the matrix elements of B u ' u, '(k). Given the m- fermion states | f m Vi), we have with 
respect to the U(N ) algebra, f m = {l m }, the antisymmetric irrep in Young tableaux notation 
and Vi are additional labels. Note that the v label used for B’s corresponds to the Young 
tableaux {2 t/ l N ~ 2l/ }. Now, Wigner-Eckart theorem for U(N) D G (here G is some subalgebra 
of U(N) giving v and ui labels) gives 


(fmVf | B v ^{k) | fmVi) = (frr 


B^k) || f m ) C^ 

JmUf 5 JmUi 


( 6 ) 


Here, (-||-1|-) is the reduced matrix element and CZZZZ is a U(N ) D G Clebsch- 

Gordan (C-G) coefficient [note that we are not making a distinction between U(N ) and SU (IV)]. 
Also, if | f r v) represents a state for r number of fermions, then | f r v) represents the corresponding 
state for r number of holes or N — r number of fermions (see [14] for details). In Young tableaux 
notation f m = {l w_m }. Definition of B v,Wv (k) and the U(N) Wigner-Racah algebra will give, 


I (fm II B u (k ) II f m ) I 2 = A U (N, m, m - k ) , 


A ti (N\ m', r ) 



N' — m! + r — n 
r 


(7) 


The A U (N, m, k) are nothing but, apart from a N and m dependent factor, a U(N) Racah 
coefficient [14]. This and the various properties of the U(N) Wigner and Racah coefficients give 
two formulas for the ensemble average of a product any two m particle matrix elements of H , 


(fmV 1 I H(k) | f m V 2 ) {f m V3 \ H(k) \ f m V 4 ) 
k 

= v 2 y A "(N,m,m-k)cy - r -^C v f v , 

n fmV 1 , f m V 2 fmV 3 , / m «4 

V=0\u>v 


( 8 ) 


and also 


{fmV 1 I H(k) I f m V2) {fmV 3 | H(k) \ f m V 4 ) 
m—k 

= v 2 y A u (N,m, k) c v, “ v 

* ^ fmV 1 5 JmV 4 JmV 3 , fm'^2 


(9) 


Eq. (9) follows by applying a Racah transform to the product of the two C-G coefficients 
appearing in Eq. (8). Let us mention two important properties of the U(N) C-G coefficients 
that are quite useful, 


V cy — 

' ^ JmVi , JmVi 
Vi 



C T T ~— 

JmVi , JmVj 



( 10 ) 


From now on we will use the symbol f m only in the C-G coefficients, Racah coefficients and 
the reduced matrix elements. However, for the matrix elements of an operator we will use m. 
implying totally antisymmetric state for fermions. An important by-product of Eqs. (9) and 
(10) is 

y {mvi I H(k) I mvj) (mvj \ H(k) \ mvk) = ([H(k)] 2 ) m 5 Vi , Vk (11) 

v i 

and we will use this in Section 3. 

Starting with Eq. (4) and using Eqs. (5), (9) and (10) will immediately give the formula, 


(MT = f‘ v ) 

\m J 


y^ (mvi | [ H(k )] 2 | mvi) 


Vl A°(N,m,k) . 


( 12 ) 




Similarly, for (H^) m first the ensemble average is decomposed into 3 terms as, 


({[H(k)] 4 )) m = 22 ( mv i I I mv i) 

Vi 

= 22 ( mv i | H(k) | mvj) (mvj \ H(k) \ mv p ) ( mv p \ H(k) \ mvi) { mv i \ H{k) \ mvi) 


Vi,Vj,V p ,Vl 


+{mvi | H(k ) | mvj) (mvi \H(k) \ mvi) (mvj \ H(k ) | mv p ) (mv p \ H(k ) | mvi) 


+ (mvi | H(k) | mvj) (mv p \ H(k ) | mvi) (mvj \ H(k) \ mv p ) {mvi \ H(k) \ mvi) 

(13) 

Note that the trace ((i7 4 )) m = (^) It is easy to see that the first two terms simplify to 

give 2 [(H 2 ) m ] 2 and the third term is simplified by applying Eq. (8) to the first ensemble average 
and Eq. (9) to the second ensemble average. Then, the final result is 


{[H(k)]±r=2 uip) 


N 
1 m 


_l min(k,m—k) 


22 A u (N,m,k) A u (N,m,m —k) d(N : is) ; (14) 


d(N : is) = 


u=0 


N 


N 


(15) 


Now, we will derive the formulas for the moments of the transition strength densities generated 
by a transition operator O that removes fco number of particles from a m-particle system. 


3. Lower-order moments of transition strength densities: results for particle 
removal operators 

Particle removal (or addition) operators are of great interest in nuclear physics. For example one 
particle (proton or neutron) removal from a target nucleus gives information about the single 
particle levels in the target and similarly, two-particle removal gives information about pairing 
force. Let us begin with a particle removal operator O and say it removes fco number of particles 
when acting on a m fermion state. Then the general form of O is, 

0 = 22v ao A ao (k 0 ) . (16) 

ao 

Here, A aQ (ko) is a ko particle annihilation operator and ao are indices for a ko particle state. 
Note that A ao (ko) transforms as {fk 0 } = {l^ - ^ 0 } with respect to U(N) and Aa 0 {ko) transforms 
as {fk 0 }- It important to recognize that the O matrices will be rectangular matrices connecting 
m particle states to m — ko particle states. In the defining space, the matrix will be a 1 x do 
matrix with matrix elements given by V ao . Note that ao takes do values and do = (^). We will 
represent O by EGUE implying that the defining space matrix elements V ao are zero centered 
independent Gaussian random variables [also they are independent of the 1 ij(k) variables in Eq. 
(2) and therefore also independent of the W variables in Eq. (4)] with variance satisfying 

VaYp = 5 aP . (17) 

In many particle spaces the O matrix will be a d\ x d 2 matrix connecting d\ = (^) number of 
m-particle states to d 2 = ( m ^ fco ) number of (m — fco)-particle states. Using Eqs. (16) and (17), 



we have 

WW=Vo(fy , JdW T =vS( N ~ o m y (18) 

Similarly, Eq. (11) gives the relations, 


{OWHP) m = (0te>) m (. HP) m , (O^HPO) 171 = (0t0) m (HP) m ~ k 0 . (19) 

Another useful result follows by introducing complete set of states between the O' and O 

operators in Eq. (18) and applying the Wigner-Eckart theorem, 

(m || A\ko) || m — koj (m — ko || A(ko ) || m) = ^ . (20) 

Following the procedure used in Section 2, it is possible to derive formulas for the lower order 

moments of the transition strength densities generated by O defined by Eq. (16). The bivariate 
moments are defined by 

M pq = (O^HQOH p ) m . (21) 

and we will consider the moments P + Q = 2 and 4 (the P + Q = 3 moments are zero as we are 
using independent EGUE representations for O and H matrices). 

Firstly, Eqs. (19) gives, 


M 20 = ( OW) m (H2) m , Mq 2 = (OW) m (H^) m - k0 , (99) 

m 40 = (e>te>) m (H 4 ) m , Af 04 = (e>to) m (H*) m - ko . 

Now, Eq. (18) along with Eqs. (12) and (14) will give the formulas for M 2 o, AI02, -M 40 and M 04 . 
Formula for the first non-trivial moment M\ \ = (Of HOH) m is derived by introducing complete 
set of states between O t and H, H and O and O and H in the trace giving, 

Q M iiM = Q WhOHT = Y. < m > «1 I ° f I m - ko, v 2 ) (m - k 0 , v 3 \ O \ m, v 4 ) 

Vl,V2 

x (m — ko, V2 | H | m — ko, V3) (m, V4 \ H \ m, v\) . 

(23) 

Using Eq. (16) and applying Eq. (17) along with Eqs. (4) - (7) and the Wigner-Eckart theorem 
will give, 


Mn(m) = UjUj- 0 1 (£:£;) X] [A u (N,m- k 0 ,m- k 0 - k) k v {N,m,m- k)} 

u=0 

^-jfk^,°L qV,UJ v qV,(jJ v 


1/2 


X 


y „ __ _ __ 

‘ fmV\ , fm—kQ V 2 fm—kQ^3 , fmV4 fm—kQ ^2 , fm—kQ ^3 fmV 4 , fmV 1 

V1,V2,V3,V4 ;a;uj u 


Simplifying the four C-G coefficients will give finally, 

Vi,(m) = VlVl Q - 1 (££) 't z n(N,m,k 1 ,k,v) 


v =0 


Zn(N, m, ko, k, u) = ( fco ) d(N : v) A U {N, m, m — k ) A U (N, m — ko, m — ko — k) 


11/2 


(24) 


(25) 


X (_l )‘Kfm,fm-k 0 ,fk 0 )+<Kfm-k 0 ,fm-k 0 ,'') U(f m f m _ ko f m f m _ ko ] /fc Q v) . 



Here 0 is a phase factor and it is a function of the U(N) irreps. It is shown elsewhere (V.K.B. 

Kota and Manan Vyas, in preparation) that (— Jk 0 )+<t>{fm-k 0 ,fm-k 0 ,v) [/(-) will 

be positive where U (-) is a U(N) [/-coefficient. Therefore we need only U' 2 and the formula 

for this is given by [18], 


LV+10 (N-v-pj , N _ 2 1) 

\tt( f ~p~ f -p . -p \i ^ _ \ v ) \p—v) \ m—p / v ' 

\y (Jm 1 Jpi Jmi Jp j Jm—p v )\ — — 


rrycj("+ 1 ) 


(26) 


Turning to the fourth order moments, we need M 13 , M 31 and M 22 - As 0 1 7 ^ 0, here 
M 13 7 / M 31 [similarly M 40 7 / A/q 4 and M 20 7 / M 02 as seen from Eq. (22)]. Following the 
procedure used for deriving the formula for Mn(m), we have for M 31 (to) 


M 31 (m) = = 2 (fr 2 ) m Mn(m) 

+ rj rl C)- 1 X 

Vl ,^2 ,^3 ,t/4 ,V5,1/6:O',1^1 jCJj/j ,1^2 ,^i/ 2 

<^to, ui | Ai(fc 0 ) | m-k 0 , v 2 ^ (m - /c 0 , ^3 I A a (k 0 ) \ m, v 4 ) 


(27) 


x (m — feo, U2 | B Ul ’ Uv i (k) | m — ko, U3) (to, v$ \ B Ul ’ bJv 1 (fc) | to, vq) 

x (m, u4 | B U 2 ,UJv ^{k) | m, u5) (to, u6 | B V 2 ,u "^{k) \ m,v 1) . 

Now, applying the Wigner-Eckart theorem, using the results in Section 2 and simplifying the 
resulting C-G coefficients will give, 

M 31 (to) = (O^HOH 3 ) m = 2 (ET 2 ) m Mu (to) 

min(k,m—k) (28) 

+vgvg (2)-‘ Cd:) X V(JV,m,t) ZnWm.ko.k,*) . 

v=0 

The function Z\\ is defined in Eq. (25). Following the same procedure as above, the formula 
for Mi 3 is, 


Mi 3 (to) = {Oni 3 OH) m = 2 (H 2 ) m ~ ko Mn( to) 

min(k,m—ko—k) 

+V&V 2 E A u (N,m-k 0 ,k ) Zii^to^o,^) • 

*z=0 


(29) 


Derivation of the formula for M 2 2 is more involved. Leaving details to a long paper under 



preparation, the final result (with p a multiplicity label) is 


M 22 (m ) = (O^H 2 OH 2 ) m = (C>t Q) m {H 2 ) m (. H 2 ) m ~ k ° 

+V£ V* ( 

x J2 ( m II [B Ul (k)B U2 {k)} u:p || m) (m-ko\\ [B ul (k)B V2 (k)] v:p || m- k 0 ) 

p 

X (_lfUrnJ m -k 0 M)+4>U m -k 0 ,f m -k 0 ,v) Uifmfm-kofmfm-ko I fk 0 v) . 

The moments Mpq can be converted into reduced (scale free) cumulants kpQ that gives 
information about the shape of the bivariate transition strength density. For our purpose the 
first non-trivial cumulants are the fourth order cumulants and they are given by, 

&40 = M40 — 3 , fco4 = Ato4 — 3 , 

kn = P31 - 3 £ , fci 3 = /ii3 — 3 ^ , 

k 22 = p 22 — 2 £ 2 — 1 ; (31) 

( r~ iP/2r~ -IQ/2'1- 1 — ~ 

MPQ = <! M 2 q M 02 > M P q and Mpq = Mpq/M 00 . 

The kpQ, P + Q = 4 follow from Eqs. (12), (14), (18), (22),(25), (28), (29) and (30). Numerical 
results for some typical values of ( N,m,k,ko ) are shown in Table 1. These results show that 
in general |Mpq| ^ 0.3 indicating that the bivariate strength density will be close to a bivariate 
Gaussian. For further confirming this result, we will derive asymptotic results for kpQ. 

4. Asymptotic formulas for bivariate moments and approach to bivariate Gaussian 
form 

Here we will consider the asymptotic limit defined by N —> oo with m , k and k$ fixed and k, 
ko « m. Note that in the dilute limit (or true asymptotic limit) we also have m —>• oo and 
m/N —> 0 with k and fco fixed. Using the formulas given in Sections 2 and 3 first we can show 
that in the asymptotic limit (asymp), 

1 Icy 

C)- 1 C- k k ° 0 ) z u ( N > m > fc o, k, k) = [(3 d(N : k) A k (N, m, m - k ) A k (N, m - k 0 , m - k 0 - A;)] 

X \U(f m ,J^r 0 J m ,frn-k 0 -,fk 0 ,k)\ aS n lP (™) (*) (V) . 

PWT = A°(N,m,k) as ^ ip (-) (£) , 

(m) _1 A fe (iV, m,m — k) A k (N,m, k) d(N : fc) aS ^ P (™) ( m “ fc ) (?)' • 



( 32 ) 



Table 1. Bivariate correlation coefficient (£) and fourth order bivariate cumulants k rs with 
r > s and r + s = 4 for various values of number of sp states (IV), number of fermions (m), 
Hamiltonian body rank (k) and the rank (/cq) of the particle removal transition operator. Results 
are obtained using the formulas given in Section 3. Note that for the A /22 that is needed for k. 22 , 
we have used Eq. (30) with the third term replaced by the corresponding asymptotic formula 
given by Eq. (37) as a formula for the reduced matrix elements in Eq. (30) is not available. 


N 

m 

k 

ko 

£ 

ki 10 

&04 

ki\ 

k\‘i 

&22 

20 

10 

2 

1 

0.82 

-0.54 

-0.55 

-0.44 

-0.45 

- 0.21 

30 

10 

2 

1 

0.85 

-0.48 

-0.50 

-0.41 

-0.43 

-0.26 

60 

10 

2 

1 

0.88 

-0.42 

-0.46 

-0.37 

-0.40 

-0.30 

80 

10 

2 

1 

0.88 

-0.41 

-0.45 

-0.36 

-0.39 

-0.31 

50 

12 

2 

1 

0.89 

-0.38 

-0.40 

-.034 

-0.36 

-0.25 


15 

2 

1 

0.91 

-0.33 

-0.35 

-0.30 

-0.31 

-0.19 


20 

2 

1 

0.92 

-0.29 

-0.29 

-0.26 

-0.27 

-0.13 


25 

2 

1 

0.92 

-0.27 

-0.27 

-0.25 

-0.25 

-0.08 

24 

8 

2 

1 

0.82 

-0.56 

-0.61 

-0.46 

-0.49 

-0.31 



2 

2 

0.66 

-0.56 

-0.67 

-0.37 

-0.43 

- 0.22 

40 

15 

2 

1 

0.90 

-0.36 

-0.37 

-0.32 

-0.33 

-0.18 



2 

2 

0.80 

-0.36 

-0.38 

-0.29 

-0.31 

- 0.12 

60 

20 

2 

1 

0.93 

-0.27 

-0.27 

-0.25 

-0.25 

-0.14 



3 

1 

0.89 

-0.51 

-0.53 

-0.46 

-0.47 

-0.30 



3 

2 

0.79 

-0.51 

-0.54 

-0.40 

-0.43 

- 0.22 



For A/ 31 , the first term in Eq. (28) is trivial and in the sum in the second term only the v = k 
term will survive in the asymptotic limit. Now, applying Eqs. (32) and (7) will give the result 



for k 3 i(m), 



Similarly ki 3 (m) is given by, 


£M = CM^ 4 oM • 


k 31 (m) 


asymp 


(m— ko—k\ (m—k\ /m\ 1/2 

1 k A k 0 M k > 

/m\ /m — ko\ 3/2 

U 0 / v k ) 


CM = CM^mM • 


(35) 


(36) 


Finally, in M 22 only the third term in Eq. (30) is complicated. This is simplified using its 
relation, valid in the asymptotic limit, to £(m) as described in [10]. Following this we have for 


&22) 


&22M 


asymp 


/m\ (m— fc \2 (m— 2 k\ fm—k\ 

9 Elm'll 2 -I- k ^ k ° ’ -I- ^ k o ^ k ’ 

T /m-ko\ /m \ 2 T (m\ (m-k 0 \ 

\ k ) \ko) W \ k ) 
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In the dilute limit with m — > 00 and m/N —>• 0 and expanding the binomials in Eqs. (33) to 
(37), it is seen that to order 1/m the cumulants k rs , r + s = 4 will be — k 2 /m (independent of 
ko) and the correlation coefficient £(m) —» 1 — (kko)/2m. Thus, the cumulants will tend to zero 
giving bivariate Gaussian form. However, as C — 1 ► 1 as m —> 00 , in practice it is necessary to 
add the k rs , r + s = 4 corrections to the bivariate Gaussian. This is same as the result seen for 
t-body transition operators before in [10, 19]. 


5. Extensions to particle addition operators and beta decay type operators 

Firstly, particle addition operator is 0 + = Yl a A' a (ko) and acting on a m-particle state it 
will generate m + ko particle states. It is easy to see that the formulation in Section 3 will apply 
directly to 0 + operator by appropriately changing everywhere m — ko by m + ko giving formulas 
for £ and M rs , r + s = 4. Explicit formulas are not given here due to lack of space. 

Now we will turn to beta decay type operator and for this consider a system with m\ fermions 
in N\ sp states and m 2 fermions in N 2 sp states with H preserving (mi,m 2 )- Then, the H 
operator, assumed to be fc-body, is given by, 

H ^= E EE V aa[l3b (i,j)AUi)Mi)Al(j)A„(j), 

i+j=k a,/3& a,b£j (38) 

Vaa:pb(i,j) = (*, a: j, a \ H \ i, /3 : j, b) . 

Here we are using Greek labels a, (5, ... to denote the many particle states generated by fermions 
occupying the orbit with N\ sp states and the Roman labels a,b,... for the many particle states 
generated by the fermions occupying the orbit with N 2 sp states. For each (i,j) pair with 
i + j = k, we have a matrix V(i,j) in the /e-particle space and we assume that the V(i,j) 
matrices are represented by independent GUE’s with their matrix elements being zero centered 
with variance, 


Vaa'.f}b(ii j) Va’a 1 :f3'b' (i 1 j ) — l^ff(bj) &ii'$ap'&cib'$fla'$ba' ■ 


(39) 



It is important to note that the embedding algebra for the EGUE generated by the action of 
the H{k ) operator on \mi,v a : rri 2 ,v a ) states is the direct sum algebra U(Ni) © U(N 2 ). Thus 
we have EGUE(A:)-[17(iV 1 ) © U(N 2 )] ensemble. A beta decay type transition operator is given 
by 

o = ^ OaaA^ko) A a (k 0 ) ; O aa = (k 0 , a\0\ k 0 , a) . (40) 

a,a 

Note that for beta decay fco = 1 and for double beta decay ko = 2 in Eq. (40). To 
proceed further, we assume a GUE representation for the O matrix in the defining space giving 

O'a^Opfi = Vq SapSab- Note that in general the O matrix is a rectangular matrix. Now, the 
ensemble averaged bivariate moments of the transition strength density are MpQ{rri \, m 2 ) = 

(P^HQOH p )^ mi,m2 \ Note that O takes (mi,m 2 ) to (mi + ko,rri 2 — ko). Formulas for the 
moments will follow by applying the formulation for particle removal operator (given in Section 
3) in m 2 space and for particle addition operator in mi space with appropriate summations over 
different parts of H . Using this, formulas are derived for Mpq with P + Q = 2 and 4 and there 
will be reported elsewhere. 

6. Conclusions and future outlook 

In this paper, we have presented exact (finite N) results for the moments of the transition 
strength densities generated by particle removal operators, using U(N) Wigner-Racah algebra 
for EGUE random matrix ensembles for spinless fermion systems. In particular, formulas for 
the moments up to fourth order are derived in detail for the Hamiltonian a EGUE(/t) and a 
k 0 number of particles removal transition operator with its structure coefficients in the defining 
spaces [V QQ in Eq. (16)] assumed to be independent Gaussian variables. Numerical results on 
one hand and the asymptotic results derived from the exact results on the other, showed that the 
fourth order cumulants approach zero in the dilute limit implying that the strength densities 
approach bivariate Gaussian form. As discussed briefly in Section 5, the formulation given 
in Sections 3 extends to transition operators that are particle addition operators and also to 
beta decay and neutrinoless double beta decay type operators. Results of the present work, the 
results reported in the Ghent meeting in June 2014 [10] where the transition operator is a t-body 
operator represented by EGUE(t) and the results (briefly reported in [11] and to be reported 
in detail in a long paper in preparation) for beta decay and double beta decay type operators 
establish clearly that the form for the bivariate transition strength densities for isolated finite 
fermion systems will be generically a bivariate Gaussian. Therefore, the bivariate Gaussian form 
with some corrections can be used in practical applications in calculating transition strengths 
in complex systems just as the corresponding results for level densities are being applied for 
calculating nuclear level densities by the Michigan group [20, 21]. Finally, further extensions of 
the present work to EGUE with U (fl) x SU (r) embedding discussed in [22] and also to boson 
systems will be important and they will have applications in mesoscopic systems [here r = 2 
with fermions is important] and Bose gases [here r = 1 and r = 3 with bosons will be important]. 
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